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Abstract. Measurements of anisotropic flow in heavy-ion coUisions provide insight into the 
early stage of the system's evolution. This proceedings presents directed and elliptic flow for 
Au+Au collisions at 39, 11.5 and 7.7 GeV, and for Cu-I-Cu at 22.4 GeV, measured in the STAR 
Experiment at RHIC. Differential measurements of directed and elliptic flow of charged particles 
as a function of centrality, transverse momentum and pseudorapidity are discussed. 



1. Introduction 

The study of collective flow in relativistic nuclear collisions has potential to offer insights into 
the equation of state of the produced matter [U [2]. Anisotropic flow is characterized by the 
Fourier coefficients [3j 

Vn = {cos n{(p - ^^Ir)), (1) 

where (j) denotes the azimuthal angle of an outgoing particle, is the orientation of the 
reaction plane and n denotes the harmonic. The reaction plane is defined by the beam axis and 
the impact parameter [3j. 

Directed flow, vi, is the first harmonic coefficient of the above Fourier expansion of the final 
momentum-space azimuthal anisotropy, and it reflects the collective sidewards motion of the 
particles in the final state. In regions that are closer to beam rapidity (y) than to midrapidity, 
directed flow is imparted very early, at a pre-equilibrium stage of the collision [H [5] , and thus 
it probes the onset of bulk collective behavior. Both hydrodynamic and nuclear transport 
models [6l [7] indicate that directed flow is a sensitive signature for phenomena related to a 
possible phase transition, especially in the general region of beam energy under investigation 
here [8]. In particular, the shape of vi{y) in the midrapidity region is of special interest because 
it has been argued that directed flow may exhibit flatness at midrapidity due to a strong, 
tilted expansion of the source. Such tilted expansion gives rise to anti-flow or a 3rd flow [9] 
component. The anti-flow is perpendicular to the source surface, and is in the opposite direction 
to the repulsive "bounce off' motion of nucleons. If the tilted expansion is strong enough, 
it can cancel and reverse the motion in the bounce-off direction and results in a negative 
vi{y) slope at midrapidity, potentially producing a wiggle-like structure in vi{y). A wiggle 
for baryons is a possible signature of a phase transition between hadronic matter and Quark 
Gluon Plasma (QGP), although QGP is not the only possible explanation [10l|TTl[9]. If strong 
but incomplete baryon stopping is assumed together with strong space-momentum correlations 



caused by transverse radial expansion, then a wiggle structure might be explained even in a 
hadronic system [12j. 

Elliptic flow, f2j is the second harmonic coefficient of the Fourier expansion above. The initial- 
state spatial eccentricity of the participant zone drives the process whereby the interactions 
produce an anisotropic distribution of momenta relative to the reaction plane. The elliptic 
momentum anisotropy saturates quite early in the collision evolution, although a little later 
than when directed flow is imparted [H [3]. Elliptic flow can provide information about the 
pressure gradients in a hydrodynamic description, and about the effective degrees of freedom, 
the extent of thermalization, and the equation of state of the matter created at early times. 
Studying the dependence of elliptic flow on system size, number of constituent quarks, and 
transverse momentum/mass, are crucial to the understanding of the properties of the produced 
matter [T]. 

2. Methods and Analysis 

In this proceedings, we report vi and V2 measurements by the STAR experiment from ^snn 
= 39, 11.5 and 7.7 GeV Au + Au and y^s^N = 22.4 GeV Cu + Cu collisions. Data were 
taken from Run 10 (2010) and Run 5 (2005). The STAR Time Projection Chamber (TPC) 
[16 was used as the main detector for charged particle tracking and second order event plane 
determination for elliptic flow analysis. The centrality was determined by the number of tracks 
from the pseudorapidity region \ri\ < 0.5. Forward Time Projection Chambers (FTPCs) were 
also used for charged particle tracking at forward rapidities. Two Beam Beam Counters covering 
3.3 < \r]\ < 5.0 were used to reconstruct the first-order event plane for directed flow analysis. 
The pseudorapidity gap between BBC and TPC allows us to reduce some of the non-flow effects. 
Elliptic flow measurements with the standard event plane method are based on the first-order 
event plane from the BBC as well as the second-order event plane based on tracks in the main 
Time Projection Chamber [16j of STAR. We also present elliptic flow based on two- and four- 
particle correlations . We analyzed minimum bias events with the primary collision vertex 
position along the beam direction (Vz) within 30,50 and 70 cm of the center of the detector 
are selected for this analysis for Au-|- Au collisions at 39, 11.5 and 7.7 GeV respectably and 30 
cm for Cu-|-Cu collisions at 22.4 GeV . In order to reject events which involve interactions with 



the beam pipe and also to minimize beam-gas interactions, the event vertex radius {^JVx + 
where Vx and Vy are the vertex positions along the x and y directions, respectively) is required 
to be less than 2 cm. For these analyses, tracks which have transverse momenta pt > 0.2 
GeV /c, pass within 3 cm of the primary vertex, have at least 15 space points in the main 
TPC acceptance {\r]\ < 1.0) or 5 space points in the case of tracks in the FTPC acceptance 
2.5 < |t/| < 4.0). Also, the ratio of the number of actual space points to the maximum possible 
number of space points for that track's trajectory was required to be greater than 0.52, which 
prevents split tracks from being counted as two tracks. The differential measurement of directed 
and elliptic flow is presented as a function of pseudorapidity (r/), transverse momentum and 
centrality. 

Figure [U left panel, shows the first-order event plane resolution as a function of centrality for 
Au -|- Au collisions at ^/sNN = 39, 11.5 and 7.7 GeV and for Cu -|- Cu collisions at ^/snn = 22.4 
GeV based on the event plane reconstructed using BBC information. The right panel shows the 
second-order event plane resolution reconstructed using charged tracks from the TPC. The errors 
shown are statistical. Generally, the first-order event plane is used for directed flow measurement 
and the second-order event plane is used for elliptic flow measurement in the standard method of 
flow analysis. STAR is well suited for these measurements, with very good event plane resolution 
down to at least 7.7 GeV. 
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Figure 1. The left panel shows the first-order event plane resolution based on information from 
the STAR BBC as a function of centrality for Au + Au collisions at ^/sN^ = 39, 11.5 and 7.7 
GeV and for Cu + Cu collisions at ^/snn = 22.4 GeV. The right panel shows the second-order 
event plane resolution based on the TPC as a function of centrality for Au -|- Au collisions at 
y^ijviv = 39, 11.5 and 7.7 GeV and for Cu -|- Cu collisions at ^/snn = 22.4 GeV. The errors 
shown are statistical. 



3. Results and Discussion 

3.1. Directed flow 
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Figure 2. Charged hadron vi vs. rj for 0-60% centrality Cu-|-Cu collisions at ^/sJfN = 22.4 
GeV obtained from the standard method based on the BBC event plane. The errors shown are 
statistical. Results are compared to vi (solid triangles) from 0-40% centrality Au-|-Au collisions 
at ^/sj^ = 19.6 GeV [2T] from the PHOBOS cohaboration. 



Figure El shows charged hadron uijBBC} in Cu-|-Cu collisions for 0-60% centrality at 



■s/snn = 22.4 GeV, compared to 0-40% centrality Au+Au collisions at y^J/vTv = 19-6 GeV. 
It is seen that the latter (PHOBOS) results are quite similar, notwithstanding the difference in 
system size, and the fact that the centrality range and beam energy are not quite the same. At 
200 GeV and 62.4 GeV, for 30-60% central collisions, we have previously reported that directed 
flow does not vary within errors between Au+Au and Cu+Cu [22]. The new finding reported 
here suggests that this behavior extends to lower energies. Figure [3] shows charged hadron vi as 
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Figure 3. Charged hadron vi as a function of pseudorapidity r] for central 0-10 %, mid central 
10-40% and peripheral 40-80 % Au+Au coUisions at = 39, 11.5 and 7.7 GeV. The errors 

shown are statistical. 



a function of pseudorapidity rj for central 0-10 %, mid central 10-40% and peripheral 40-80 % 
Au+Au collisions at ^/sNN = 39, 11.5 and 7.7 GeV. The errors shown are statistical. We observe 
centrality dependence of the directed flow signal specially prominent at mid rapidity( < 1.0) 
than at forward rapidities. Figure U] shows charged hadron vi as a function of transverse 
momentum for central 0-10 %, mid central 10-40% and peripheral 40-80 % Au+Au collisions 
at -y/s/vjv = 39, 11.5 and 7.7 GeV. The errors shown are statistical. For this plot, particles from 
the main TPC are used with \r]\ < 1.0 . We observe the centrality dependence of directed flow 
as a function oipx ■ We also observe that the directed flow signal crosses zero at px > l.OGeV/c 
from negative to positive , which was previously reported at higher RHIC energies [22] 

Figure [5] (left panel) shows charged hadron vi as a function of rj scaled by the corresponding 
ybeam for Au+Au Collisions at 7.7, 11.5, 39, 62.4 and 200 GeV for 30-60 % central collisions. 
The new results reported here are the charged hadron wijBBC} in Au+Au collisions for 30- 
60% centrality at ^/sNN = 39,11.5 and 7.7 GeV. On the right, we plot charged hadron vi as a 
function of rj scaled by the corresponding ?/beam 

for Cu+Cu colhsions at 22.4, 62.4 and 200 GeV. 
The results for 62.4 and 200 GeV are for 30-60% centrality previously reported by STAR [22] . 
The new results reported here are the charged hadron t>i{BBC} in Cu+Cu collisions for 0-60% 
centrality at y^i^wv = 22.4 GeV. We observe that t'i(f?/ybeam) shows a beam energy scaling 
behavior, though not perfect, that has already been observed at higher RHIC energies. 

Figure [6] (left) shows charged hadron vi as a function of - ybeam for Au+ Au Collisions at 
7.7, 11.5, 39, 62.4 and 200 GeV for 30-60 % central collisions. On the right, we show charged 
hadron vi as a function oi rj — ybeam 

for Cu+Cu colhsions 22.4, 62.4 and 200 GeV. The results for 
62.4 and 200 GeV are for 30-60% centrality, previously reported by STAR [22]. In this frame. 
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Figure 4. Charged hadron vi as a function of transverse momentum px for central 0-10 %, mid 
central 10-40% and peripheral 40-80 % Au+Au collisions at ./s^ = 39, 11.5 and 7.7 GeV. The 
errors shown are statistical. For this plot, particles from the main TPC are used with |r/| < 1.0 
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Figure 5. The left panel shows charged hadron vi as a function of r/ scaled by the respective 
T/beam for Au+ Au Collisions at 7.7, 11.5, 39, 62.4 and 200 GeV for 30-60 % central collisions. 
The right panel shows charged hadron vi as a function of t] scaled by the corresponding ybeam for 
Cu+Cu cohisions at 22.4, 62.4 and 200 GeV. The results for 62.4 and 200 GeV are for 30-60% 
centrality, previously reported by STAR [22] . 



zero on the horizontal axis corresponds to beam rapidity for each of the these beam energies. 
The data support the limiting fragmentation hypothesis [22] in the region —2.0 < ybcam < — 1 
for both system Au+ Au and Cu+ Cu presented here. 
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Figure 6. The left panel shows charged hadron vi as a function of rj — ybcam for Au+Au 
collisions at 7.7, 11.5, 39, 62.4 and 200 GeV for 30-60 % centrality. The right panel gives 
charged hadron Vi as a function of TJ — ybeam 

for Cu+Cu colhsions at 22.4, 62.4 and 200 GeV. 
The results for 62.4 and 200 GeV are for 30-60% centrality, previously reported by STAR [22j. 



3.2. Elliptic flow 
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Figure 7. The left panel shows elliptic flow V2 as a function of centrality for charged hadrons 
from Au+Au collisions at ^/sJlN= 7.7, 11.5, 39, 62.4 and 200 GeV. The right panel shows 
elliptic flow V2 as a function of centrality for Cu+Cu collisions at -^/ijvA''=22.4, 62.4 and 200 
GeV. These results are based on the two-particle direct cumulant method with \r]\ < 1.0 and 
0.2 <PT < 2.0GeV. 



Figure [7] (left) shows elliptic flow V2 as a function of centrality for charged hadrons from 
Au+Au collisions at yJSNN= 7.7, 11.5, 39, 62.4 and 200 GeV. The right panel presents elhptic 
flow V2 as a function of centrality for Cu+Cu collisions at y^i7V7V=22.4, 62.4 and 200 GeV. These 
results are based on the two-particle direct cumulant method with |r/| < 1.0 and 0.2 < pT < 2.0 
GeV. We observe that V2 as a function of centrality increases with beam energy. 
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Figure 8. The left panel shows V2 as a function of rj for Au+Au collisions at 7.7, 11.5, 39, 
62.4 and 200 GeV for 0-40 % central collisions. On the right, we show the elliptic flow versus r/ 
for charged hadrons from 22.4 GeV Cu+Cu collisions at 0-60% centrality measured with the rj- 
subevent method with an ij gap of 0.3 units. Results are compared with data from the PHOBOS 
collaboration at 22.4, 62.4 and 200 GeV with 0-40 % centrality. 



Figure [5] left panel shows V2 as a function of r] for Au+Au collisions at 7.7, 11.5, 39, 62.4 and 
200 GeV for 0-40 % centrality. On the right, we show elliptic flow versus rj for charged hadrons 
from Cu+Cu collisions at 0-60% centrality at ^/sJfN = 22.4 GeV compared with the results 
from the PHOBOS collaboration at 22.4, 62.4 and 200 GeV with 0-40 % centrality. The new 
22.4 GeV Cu+Cu results reported here are based on the r/-subevent method with an rj gap of 0.3 
units. We observe that f 2 as a function of rj increases with beam energy. This energy dependence 
of pt integrated elliptic flow as a function of centrality and rj is driven by the increase in mean 
Pt with energy. 

Figure [9] left panel shows V2 as a function of px for Au+Au collisions at 7.7, 11.5 and 39 
GeV for 0-80 % centrality. The results reported here are based on the event plane method 
using full event plane from TPC and are consistent with the result obtained with two particle 
cumulant methods. On the right, we show elliptic flow versus pT for charged hadrons from 
Cu+Cu collisions at 0-60% centrality at ^/snn = 22.4 GeV, compared with the previously 
published results at 200 GeV. Again, the new 22.4 GeV Cu+Cu results reported here are based 
on the 7/-subevent method with an r] gap of 0.3 units. The error bars include only statistical 
uncertainties. We observe V2 as a function of obtained using two particle correlation as well as 
event plane method for Au+Au collisions shows the energy dependence at lower beam energies 
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Figure 9. The left panel shows V2 as a function of px for Au+Au collisions at 7.7, 11.5 and 
39 GeV for 0-80 % centrality. The right panel shows elliptic flow versus pT for charged hadrons 
from Cu+Cu collisions at 0-60% centrality at y/s^N = 22 A GeV measured with the 77-subevent 
method with an rj gap of 0.3 units. Results are compared with previously published data at 200 
GeV. 



upto 39 GeV. However the result obtained using four particle cumulant method [29] shows very 
small energy dependence. 

Figure [10] presents elliptic flow W2{4} for charged hadrons as a function of transverse 
momentum for various beam energies from 7.7 GeV to 2.76 TeV. The results from 7.7 to 200 
GeV are from Au+Au collisions at RHIC and those for 2.76 TeV are from Pb+Pb collisions at 
the LHC [30] . The red line shows a fit to the results from Au+Au collisions at ^/sNN= 200 GeV. 
The bottom panel shows f2{4} divided by the fitted line, plotted vs. px for all energies. The 
results are shown for 10-20 % central collisions. It is observed that f2{4}(pT) for all energies 
shows similar values beyond pT ~ 500 MeV/c. The agreement is within ~ 10 % of the 200 GeV 
data if we consider only energies of 39 GeV and above. The statistical errors at the energies 
of 11.5 and 7.7 GeV are larger. This saturation is remarkable considering the very wide energy 
range involved. 

4. Summary 

In this proceedings, we present STAR results for directed and elliptic flow as a function of 
transverse momentum, pseudorapidity and centrality for Au+Au collisions at 39, 11.5 and 7.7 
GeV and for Cu+Cu collisions at 22.4 GeV. Our findings extend to 22.4 GeV the observation 
that directed flow is independent of system size (comparing Cu+Cu with Au+Au) at 62.4 and 
200 GeV. Our findings also demonstrate that vi(r//2/beam) shows a beam energy scaling behavior, 
though not perfect, that has already been established at higher RHIC energies. When viewed 
in the projectile frame of reference, our results are consistent with the limiting fragmentation 
hypothesis. Elliptic fiow as a function of centrality and pseudorapidity increases with beam 
energy in both systems studied, i.e., Au+Au and Cu+Cu. The centrality and pseudorapidity 
dependence are similar to that observed at higher RHIC beam energies. The measured elliptic 
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Figure 10. Elliptic flow U2{4} for charged hadrons as a function of transverse momentum 
for various beam energies from 7.7 GeV to 2.76 TeV. The results for 7.7 to 200 GeV are from 
Au+Au collisions and those for 2.76 TeV are from Pb+Pb collisions [30]. The red line shows a 
fit to the results from Au+Au collisions at ^/sNN= 200 GeV. The bottom panel shows 1^2 {4} 
divided by the fitted line, plotted vs. pT for all energies. The results are shown for 10-20 % 
central collisions. 



fiow as a function of transverse momentum increases from 7.7 to 39 GeV but shows saturation 
above that energy up to 2.76 TeV. 
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